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Abstract 



CN ' The cosmology of the Randall-Sundrum scenario for a positive tension brane 

qq ' in a 5-D Universe with localized gravity has been studied previously. In the 

ON ■ radiation-dominated Universe, it was suggested that there are two solutions 

for the cosmic scale factor a(t) : the standard solution a ~ t 1 ' 2 , and a solution 

a ~ i 1 / 4 , which is incompatible with standard big bang nucleosynthesis. In 

this note, we reconsider expansion of the Universe in this scenario. We derive 

and solve a first order, linear differential equation for H 2 , the square of the 

expansion rate of the Universe, as a function of a. The differences between 

JL , our equation for H 2 and the relationship found in standard cosmology are (i) 

D \ there is a term proportional to density squared (a fact already known), which 

is small when the density is small compared to the brane tension, and (ii) 

there is a contribution which acts like a relativistic fluid. We show that this 

^ ■ second contribution is due to gravitational degrees of freedom in the bulk. 

Thus, we find that there need not be any conflict between cosmology of the 

Randall-Sundrum scenario and the standard model of cosmology. We discuss 

how reheating at the end of inflation leads to the correct relationship between 

matter density and expansion rate, H 2 — > &irGp m /?> : and the conditions that 

must be met for the expansion rate of the Universe to be close to its standard 

model value around the epoch of cosmological nucleosynthesis. 
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I. INTRODUCTION 

Recently Randall and Sundrum M presented a new static solution to the 5-D (classical) 
Einstein equations in which spacetime is flat on a 3-brane with positive tension provided 
that the bulk has an appropriate negative cosmological constant. Even if the fifth dimension 
is uncompactified, standard 4-D gravity (specifically, Newton's force law) is reproduced on 
the brane. In contrast to the compactified case ||, this follows because the near-brane 
geometry traps the massless graviton. 

To see if such a scenario is viable phenomenologically, one application to check is the 
evolution of the early Universe [|3],[|] . It was pointed out [||] that a 5-D Universe with branes 
may have a rather unconventional, and perhaps unacceptable, cosmology. The inclusion of 
matter inside the brane is suggested naturally by the brane world picture. So it is important 
to see if the standard expanding Universe can be recovered by extending the static solution 
to a time- dependent one when matter /radiation inside the brane is included. As shown in 
Ref |J, the standard matter-dominated expanding Universe is recovered for large enough 
brane tension. However, Ref |J proposed two possible behaviors for the cosmic scale factor, 
a(t), during the radiation-dominated Universe: (i) a(t) ~ t 1//4 (first found in Ref ||), or (ii) 
a(t) ~ t l l 2 (found in Ref ||). That there might be two different powerlaw solutions is not 
surprising, since a non-linear second order differential equation may have more than one 
such solution, and initial conditions will determine the precise time evolution of a(t) (which 
may not be exactly a powerlaw, but might be close to either of the candidate powerlaws 
at different times). On the other hand, to match the known observations of the expanding 
Universe, at least back to the time of electron-positron annihilation and nucleosynthesis, the 
expansion rate of the Universe should be approximately H 2 = 8irGp/3, its value in standard 
Big Bang cosmology. Actually this requirement is stronger than demanding a(t) ~ t 1 ^ 2 : 
agreement with the well-established picture of light element synthesis in the early Universe 
constrains if as a function of temperature T for T ~ 100 keV - a few MeV. 

In this note, we clarify the situation with respect to the cosmology of the Randall- 
Sundrum scenario. First, we elucidate why the evolution equation that one obtains for 
the scale factor [Eq. ( P^ ) below] is second order in time rather than first order in time 
as in standard cosmology. The reason is that a piece of the gravitational dynamics in 
the bulk is coupled to the brane dynamics, and thus there is an extra free parameter in 
the cosmological equations, describing the amount of effective 4-D energy density due to 
the 5-D gravitational degrees of freedom. The most general bulk metric compatible with 
homogeneity and isotropy is a black hole anti-deSitter metric ||, and the free parameter is 
determined by mass parameter of this bulk solution. This conclusion has also been reached 
by Kraus ||. 

Second, we derive and solve a first-order differential equation for H 2 as a function of 
log a. This equation shows that the effect of bulk gravitational degrees of freedom on the 
4-D dynamics is to add a new term that decays with expansion just like a relativistic fluid 
(but which may be positive or negative). This new contribution to H 2 is degenerate with the 
contributions from relativistic particles such as photons and neutrinos in the early Universe. 
Thus, if significant, it would alter the relationship between H and temperature T of the 
early Universe, even though it would still allow a(t) ~ t 1 / 2 . Since it is the dependence of H 
on T that is constrained by comparing observed light element abundances with the theory 



of cosmological nucleosynthesis 0,0' the new terms due to bulk gravitational degrees of 
freedom must be small enough to preserve the concordance between theory and observation. 
We show that if we start with an inflationary epoch [|TT1 ,|3|,|i"2f , and follow the Universe as 
it evolves into the radiation-dominated phase, it is possible (even likely) that the new term 
becomes insignificant by the time the Universe reheats at the end of inflation. Our final 
result for H 2 shows that standard cosmology can be recovered at late times in the Randall- 
Sundrum scenario in both the radiation and matter dominated eras. 

II. DERIVATION OF THE BRANE DYNAMICS 



In this section we derive the equations governing the dynamics of matter and geometry 
on the brane, using a local expansion of the metric near the brane. We consider a 5-D 
spacetime with coordinates x A = (x° , x l , x 2 , x 3 , x 4 ) = (t,^ 1 ,^ 2 ,^ 3 ,?/), and we assume that 
there is a single brane located at y = 0. We assume the following form for the 5-D metric [] 



ds 2 = ex.p{2j3(y,t)](-dt 2 + dy 2 ) + exp[2a(y, t)]8 ij dx' l dx j . 



(1) 



Thus, the metric only depends on t and y, and is flat in ordinary 3-D space (labeled by Latin 
indices which run over 1,2,3). For simplicity, we also restrict attention to the case where 
a and (3 are even functions of y. For this metric the non-zero components of the Einstein 
tensor, as shown by Binetruy et al. [H, are 
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where dots denote derivatives with respect to t and primes with respect to y. 

We assume that a and j3 are smooth functions of \y\ and of t, i.e., a(y, t) = a(\y\,t) and 
/3(y,t) = f3(\y\,t), where the functions a(£,t) and /?(£,£) are smooth in a neighborhood of 
£ = 0. Then the derivative da/dy will generically be discontinuous across the brane, as in 
Ref. [HI. We define 



1 Binetruy et al. used the metric 

ds 2 = -n 2 (y,t)dt 2 + b 2 (y,t)dy 2 + a 2 (y,t)5 ij dx i dx j , 

but did not employ the freedom to redefine the coordinates y and t to put the metric in the simpler 
form (||) used here. 
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together with a similar equation for /3 w . 

We now substitute the relation (||) into the Einstein tensor components (||) — (||), and 
insert into the 5-D Einstein equations Gab — ^ 2 Tab where the energy momentum tensor^ 
is 

T A B = diag[(-p,p,p,p,0)<%)exp(-/3) + (A 6 , A 6 , A h , A 6 , A 6 )]. (9) 

Here p and p are the density and pressure of the matter on the brane, and A b > is 
proportional to minus the bulk cosmological constant 0. The result is, first, two equations 
obtained by equating the coefficients of the S(y) in the Goo and Gij equations: 

2 

«i = -^ exp[/5(0,t)] (10) 

A= K2( Y 3p) exp[mt)]. (11) 

Second, there are the smooth pieces of the equations. Since we can restrict attention to y > 
for these pieces of the equations (by the evenness assumption), we can drop the distinction 
between a and a, and between £ = \y\ and y, so we replace terms like a^(\y\,t) [cf. Eq. (||) 
above] by a m (y,t). The resulting equations are, for y > 0, 

a 2 + a/3 - a" -2a' 2 + a' (3' = -—A b e 2(3 , (12) 

- 2d - 3d 2 - /3 + 2a" + 3a' 2 + /T = K 2 A 6 e 2/3 , (13) 

-d-2d 2 + d/3 + a /2 + a'/3' = — A b e 2f3 , (14) 



and 



/3'd + a'/3 - a' - aa' = 0. (15) 



2 The representation of the brane stress-energy tensor as a <5-function in y involves the tacit as- 
sumption that the thickness of the brane is smaller than |/3i| _1 or |ai| _1 . 

3 The negative bulk cosmological constant ensures that one solution of the bulk equations is anti- 
deSitter spacetimc M. 



A. Power Series Expansion for the Metric Near the Brane 



The equations ( 110 ) - (|T5|) define the coupled dynamics of the brane and bulk degrees of 
freedom. In this section we derive a description of the dynamics of the brane by itself, by 
using a power series expansion for the metric near the brane. 

We assume that, near y = 0, 

a(y,t) = a (t) +ati(t)\y\ + -a 2 (t)y 2 H 

/3(y, t) = fat) + M)\y\ + \(3 2 {t)y 2 + ■■■■ (16) 

This assumption is compatible with our definitions (0) and (0) of a,\ and p\. If we insert 
this expansion into Einstein's equations we obtain once again the relations (fLOl) and (JTTJ) 
on the brane. Equating powers of the smooth pieces ( |T2| ) - (|T5j) of the equations of motion 
yields, first, from the piece of the G04 equation ( [Llf) that is ex y°, 

p + 3d (p + p) =0, (17) 

and five others that follow from the remaining terms, 

•2 , a ■ o 2 , o K 2 A b exp{2f3 ) 
a + p ao - 2a x + aipi — a 2 = - (18) 



- 2a - 3dg - fa + 3aj + 2a 2 + (3 2 = K 2 A b exp(2(3 ) (19) 

- a 2 + oc 2 ($o - do) + /?2«o + A«i + $iOii - di«i = (20) 

o-2 , • a . / .on K 2 A&exp(2/5 ) , . 

-a -2a + a /5 + ai(ai + p\) = (21) 



^o , a \ i a - a • ■ i • a i « • 2k 2 A 6 /5i exp(2/? ) ,„„, 

a 2 (2ai + pi) + aip 2 — «i — 4a ai + a Pi + Po«i = ^ • {^^) 

Equation ( I7|) is simply the usual conservation of energy equation, and its derivation from 
the G04 component of the 5-D Einstein equations was already given by Binetruy et al. || . 
Using Eq. (|TJ) in Eq. (|21|) we find 

- do - 2al + a p = (*£* + /tV(P 36 +3p) ) exp(2/3 ); (23) 

defining a new time variable t by dt = exp(f3 )dt, we obtain the equation already found by 
Binetruy et al. and Csaki et al. ||, 

d 2 a _ /da \ 2 _ n 2 A b n 4 p(p + 3p) 



dt 2 V dt J 3 36 



Note that i is just the conventional cosmo logical proper time, that is, the proper time 
as measured by co moving observers on the brane. If we assume that p = o + p m and 
p = —o + p m , where a represents the contribution from the brane tension, and p m and p m 
are the density and pressure due to the matter, then this equation reduces to 

d 2 a /da \ 2 _ K 2 A b k 4 o 2 K*a(3p m - p m ) n 4 p m (p m + 3p m ) 
di 2 V dt ) ~ 3 18 36 36 ' [ ' 

We shall find exact solutions of Eq. ( pop below. 

To complete the solution, note that Eqs. (|T8|) , fllTf ) and (|^) are algebraic equations for 
a 2 and /3 2 . Equation ( |T8"D implies that 

2 , • a n 2 , _ . a , ^ 2 A fe exp(2/3 ) 
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= d^ + d /5o + [-^ g 6 7 j exp(2/3 ), (26) 

where Eq. ( |TT| ) was used to get the second expression. Using this result for a 2 in Eq. flI9p 
gives 

/9 2 = /t2Afee 3 P(2/3 ° ) + «o + 2«o + A, + «i - 2«i/?i = -3dg + /3 + (^ - ^) exp(2/? ), 

(27) 

where we used eqs. ([TT]) and (21|) in the last line. It is easy to check that once we have 
determined a 2 and /3 2 , the remaining equations ( PU[ ) and ( ^2] ) are both satisfied identically. 
Note that the function /3o(t) is not determined from this power series expansion near the 
brane. Thus, it is not specified by physics on the brane itself. It is determined by the 
dynamics in the bulk, but is gauge dependent. 

III. AN EQUATION FOR H 2 

We now turn to the task of solving the second order dynamical equation (p5|) for ao(t). 
The easiest way to do this is to rewrite it as a first order equation for H 2 = (dao/dt) 2 as a 
function of «o- This rewriting may be accomplished by noting that 

d 2 a d(dao/dt) dao TjdH ldH 2 
di 2 da di da 2 da ' 

then Eq. (p5|) becomes 

^dH 2 2 exp(-4a ) d[H 2 exp(4a )} 
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(30) 



Note that Eq. (|30"D represents a linear, first-order equation for H 2 , if we regard the right 
hand side as a source term that is (implicitly) a function of o,q. Implicit in our change of 
variables from t to ao is the restriction to phases of the evolution of the Universe in which 
the two variables are related to one another monotonically. For oscillating Universes, or 
Universes in which the scale factor a(t) = exp[«o(^)] m ay have a contracting phase as well 
as one or more expanding ones, we can derive Eq. (|30"D for each of these phases separately. 
To solve Eq. ([30]) generally, let us first rewrite the equation of energy conservation, Eq. 
(0), as 
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this equation implies 
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Substituting this result for p m into equation Eq. (|30|) yields 
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we can easily read off the solution 
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K exp(— 4a ) 



(34) 



where K is a constant of integration, and may be positive or negative. In the case of the 
Randall- Sundrum static solution |l|, we have n 2 a 2 = 6A fc , kV = 487rG, and p m = K = 0, 
giving H = and a(t) = exp[a (t)] = constant. An equation similar to Eq. ([34]) has been 
derived by Kraus using a different approach [Eq. (27) of Ref. 0]. 

Just as in standard cosmology based on 4-D general relativity, the evolution of the Uni- 
verse can be reduced to the solution of an equation for H 2 and an equation for energy 
conservation. However, the equation for H 2 has a different structure than in standard cos- 
mology for three reasons. First, there are two terms that result from the brane tension and 
the negative cosmological constant in the bulk; in the Randall-Sundrum scenario, these terms 
may be chosen to cancel exactly. Second, in addition to the normal term proportional to the 
matter density p m , there is a term that is proportional to p 2 m . This high-density "correction" 
term, which is reminiscent of the small-scale deviation from Newton's law found in this sce- 
nario, becomes unimportant once p m <^C a. Third, there is the term i^exp(— 4a ) that arises 
purely from initial conditions. This is a qualitatively new feature of the Randall-Sundrum 
scenario. In cosmology based on 4-D general relativity, H 2 is completely determined by the 
energy density of the Universe (presuming a spatially flat model, as was done above). But 



in the reduction from five dimensions to four dimensions done here, we find that H 2 can be 
specified freely at some initial time. The additional term that results decays exactly as any 
relativistic matter density would (although it need not be positive). 

The physical interpretation of the i^exp(— 4«o) is as follows. The most general solution 
of the bulk equations of motion QT2J) — (|T5| ) is the black hole solution || 

ds 2 = -w^dt 2 + w{r)- l dr 2 + r 2 8 ij dx i dx j , (35) 

where 

w(r) = -A b r 2 + ^. (36) 

6 r z 



The parameter K is determined by the mass parameter m of the solution (35); see Eq. (27) 



of Ref. |J. Positive K corresponds torn > and negative K to m < 0. Although the 



solution fl35| ) is static, the motion of the brane through the bulk need not respect the t-time 
translation symmetry, and thus there is a contribution to the effective 4-D energy density 
which varies with brane proper time 0. 

In specific cosmological scenarios, the additional term can be important or negligible, 
depending on the magnitude of K. More concretely, let us imagine that the Universe begins 
with an early, inflationary phase of expansion, during which the energy density and pressure 
are dominated by an inflaton field, and that once the inflaton approaches its potential 
minimum, its energy dissipates into relativistic particles. For example, in one model for 
inflation and reheating (e.g. [ |i~2|j ) 

2 
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where the inflaton field obeys the equation 

d 2 (j) . ozj d<f> dV((f)) 



rfp +(3„ + r,)^ + -^ = 0, (38) 

and the radiation energy density p r is determined from 

$ + «»-r,(*)'. (30, 

Here T^ is a decay rate that leads to the production of relativistic particles, primarily from 
oscillations of the inflaton field about its minimum. It is straightforward to set up and 
solve these equations along with Eq. (|25|) numerically, and we have done so. To begin the 
integration, we must specify not only initial values of 0, d<p/dt and p r , but also the starting 
value of H. Our (limited) exploration^] of numerical solutions shows that, as expected, Eq. 



4 We chose d(/)/dt = and p r = initially in all of our numerical integrations, and also focussed on 
cases where the Randall-Sundrum condition, K 2 a 2 /6At, = 1 is satisfied. A simple effective potential 



of the form V(4>) = Vo[ ex P( — e ) — ex P( — y f2 + 2 /0o)] was assumed, with e chosen to be small but 
nonzero to guard against pathologies in dV{4>)/d<j) near = 0. Although Vo/c is a parameter that 
can be chosen arbitrarily, so far we have considered Vq/ct < 1. 



([34]) is satisfied (to the accuracy of our numerical integrations) during both the inflationary- 
epoch and the radiation-dominated era that follows reheating, even though our solutions 
were based on Eq. ( p5|) (with independent variable a instead of £), not Eq. (|3^|). 

Using Eq. (|34D along with eqs. (|37|) , ( |38"D and (|39|) , we can gain insight into how inflation 
might be affected by the new features of our relationship for H 2 . We assume that the 
Randall- Sundrum condition, K 2 a 2 /QAb = 1, holds. Let us first consider the effect of choosing 
the initial value of H, which translates (with suitable definition of cto = 0) into choosing a 
value of K. To keep matters as simple as possible, let us assume that Vo, the value of V((j)) in 
its flat portion (where we presume the inflaton starts), is smaller than a, so we can neglect the 
nonlinear term in Eq. fl34|) for H 2 . Then if the initial value of K is positive and ~ 87rGVo/3, 
the starting value of H 2 in inflation theory based on 4-D general relativity, deviations in 
H 2 from its standard cosmological value will damp away exponentially as the inflaton rolls 
toward its minimum. Given enough expansion, i^exp(— 4a ) becomes negligible compared 
with 87rGVo/3 by the end of inflation, and reheating results in a radiation density that 
dominates the i^exp(— 4«o) term in Eq. (|34"D thereafter. If the initial value of H is far larger 
than in standard cosmology, the inflaton hardly rolls at all until the Universe expands enough 
that f^exp(— 4ao) becomes comparable to 87rGVo/3; the subsequent evolution is identical to 
what would transpire if we had started with H 2 ~ 87rGVo/3. This scenario can be thought 
of as including a pre- inflation, "radiation" dominated era, followed by conventional inflation 
and reheating. If the initial value of H is well below its standard value initially, which 
means that K < (but not so small that H 2 < 0), then the total number of e-foldings 
during inflation is diminished, but can still be large enough that K exp(— 4«o) becomes 
unimportant by the time reheating occurs. Excluding negative H 2 amounts to limiting 
— K = \K\ to values smaller than k 4 Vo(ct + Vo/2)/18 initially (at some small but nonzero 
starting value of a(t); spatially flat models with K < have no Big Bang). This might not 
be unnatural, as we might expect Vo ~ a ~ \K\Jk a o. 

The nonlinear term in Eq. (Kl) can be important if Vo ^ o. In that case, the Universe 
expands faster than in standard cosmology, but H is still time-independent while <fi is rolling 
down the flat part of V(4>), so much of the standard picture of inflation can be carried over 
intact. In particular, as discussed in the previous paragraph, it is possible that the Universe 
expands enough during inflation to render i^exp(— 4ao) negligible before reheating. The 
number of e-foldings of a(t) between the time when the presently observable portion of 
the Universe crossed the horizon and the end of inflation is Nh ~ HTq/V Hq, where 
To and H$ are the present temperature and Hubble parameter. Since Nh oc H, larger H 
means more e-foldings, and, generically, larger primordial density fluctuations, which could 
be problematic. When reheating is complete, it is still possible that p r > o for awhile, until 
expansion can reverse the inequality. During the phase in which H 2 « /t 4 p^/36, the scale 
factor behaves as a{t) = exp[«o(£)] ~ £ 1//4 as found in Refs. 0,0, but once H 2 — > K 4 crp m /18, 
the scale factor behaves as a{t) ~ i 1 / 2 . As long as a is large compared with ~(MeV) 4 , the 
nonlinear term in Eq. ([34]) has no effect on cosmological nucleosynthesis. 

In the scenario outlined above, inflation plays a central role in guaranteeing that the evo- 
lution of the Universe tends toward the standard cosmological model by the nucleosynthesis 
epoch, at T ~ 100 keV - a few MeV. It is not sufficient to require that a{t) —> t 1 ' 2 at late 
times to obtain agreement with standard model predictions for the synthesis of the light 
elements, since this solution is obtained for K ^ 0, even if i^exp(— 4ao) > k 4 <tVo/18, and 



H 2 m K exp(— 4«o)- The outcome of cosmological nucleosynthesis depends on the compe- 
tition between H and various particle reaction rates, such as those that create and destroy 
protons and neutrons; these reaction rates can be computed as functions of T and the baryon 
density of the Universe [f7 ,T0f1 ■ The comparison of the observed light element abundances 



with the theoretical calculations constrains the ratio of p m to its standard model value at 
T ~ 100 keV - 1 MeV, which includes contributions from photons, three neutrino flavors, 
and, at T > m e , electrons and positrons to be very close to one [|TUJ. Thus, there is little 
room for additional relativistic species, such as K exp(— 4a ) for positive K. If K > 0, it is 
crucial that K/(a(t)T) A become negligible well before nucleosynthesis occurs, as it is (vir- 
tually) unaffected by expansion alone. If i^exp(— 4cto) < 0, it would have to be balanced 
nearly exactly by additional particles that are relativistic during the nucleosynthesis era. 
While one cannot rule out such a perverse conspiracy a priori, it would require fine tuning 
that seems rather unlikely, so even for negative K it seems that, generically, K / (a(t)T) 4 
must be negligible at T ~ 100 keV - 1 MeV. Inflation and reheating provide a natural 
mechanism for decreasing K/(a{t)TY to an acceptable level. Indeed, one might say that, if 
the Randall-Sundrum brane world scenario proves correct, inflation is required to preserve 
the success of the standard Big Bang cosmology in accounting for the abundances of the 
light chemical elements. 

IV. CONCLUSIONS AND REMARKS 

We have re-examined the evolution of homogeneous, isotropic cosmologies in the Randall- 
Sundrum brane world scenario, building on the earlier work of Binetruy et. al. [Bfl and Csaki 
et. al. ||. We derived the explicit form of the equations that governed the dynamics of 
both the bulk and brane degrees of freedom [Eqs. (|10D , (|TTD, ( |TTD and ([12]) — (|TJ) above]. 
We derived a generalization [Eq. ([34]) above] of the standard relation between the Hubble 
parameter H and energy density p m . This equation contains a term (an effective 4-D energy 
density) that acts like a free, relativistic species, and is the only effect of the 5-D gravitational 
degrees of freedom on the bulk dynamics. 

The cosmological equation we find for H contains extra terms compared to 4-D general 
relativity (a term proportional to density squared and the term due to 5-D gravitational 
degrees of freedom) but, generically, these become unimportant at late times. In particular, 
we analyzed how reheating at the end of inflation smoothly matches onto the standard 
behavior a{t) oc t 1 ! 2 of the radiation dominated epoch, and not the behavior a(t) oc t 1 / 4 
found in Ref. ||. We show that inflation can drive the expansion rate of the Universe 
toward its standard cosmological value, H 2 = 8irGp m /3 after reheating under a variety of 
conditions. Thus, the Randall-Sundrum scenario is compatible with standard cosmology. 



A recent paper by Shiromizu, Maeda and Sasaki |13j derives effective Einstein equations 
on the brane for a general bulk metric. Their independent analysis is compatible with our 
results; they find corrections to the Einstein equation that scale as density squared, and also 
a correction term E^ due both to bulk gravitational waves and to non-radiative bulk degrees 
of freedom. Bulk gravitational waves do not arise in our analysis because of the symmetries 
we assumed in the starting metric ([I]). On the other hand, our term i^expf— 4ao(t)] in Eq. 
(0) corresponds to the non-radiative piece of the tensor E^ u of Eq. (17) of Ref. |jl3fl . 



Finally, we remark that to solve the hierarchy problem, it was suggested [).l\ that we live 
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in a probe brane (or "TeV" brane), which is at a distance yo away from the Planck brane, 
where the graviton is trapped. In this case, the tension of the "TeV" brane is substantially 
smaller than that of the Planck brane. In the cosmological setting, depending on the details, 
the p 2 m term contribution to H 2 may no longer be negligible on the "TeV" brane, and the 
standard cosmology will be modified accordingly. It will be interesting to see if cosmological 
constraints can be consistent with this solution to the hierarchy problem. 
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